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Electric dipole moments constitute a competitive method to search for new physics, being par- 
ticularly sensitive to new CP-violating phases. Given the experimental and theoretical progress in 
this field and more generally in particle physics, the necessity for more reliable bounds than the 
ones usually employed emerges. We therefore propose an improved extraction of the electric dipole 
moment of the electron and the relevant coefficient of the electron-nucleon coupling, taking into 
account theoretical uncertainties and possible cancellations, to be used in model-dependent analy- 
ses. Specifically, we obtain at 95% C.L. \d e \ < 0.f4 x f0 -26 ecm with present data, which is very 
similar to the bound typically quoted from the YbF molecule, but obtained in a more conservative 
manner. We examine furthermore in detail the prospects for improvements and derive upper limits 
for the dipole moments of several paramagnetic systems presently under investigation, i.e. Cesium, 
Rubidium and Francium. 



I. INTRODUCTION 

Despite the tremendous success of the Standard Model 
(SM), the arguments for the necessity of an extension 
are compelling. Specifically, Sakharov's conditions [1] re- 
quire the presence of additional sources for CP violation 
with respect to the SM, given the observed baryon asym- 
metry of the universe. Assuming CPT invariance, elec- 
tric dipole moments (EDMs) are highly sensitive probes 
for new CP- violating phases. This renders them a com- 
petitive tool in the search for new physics (NP), com- 
plementary to both, direct searches at the LHC and 
Tevatron as well as indirect ones in flavour-changing pro- 
cesses. 

As interface between a given theory and experiment 
typically an effective Hamiltonian is used. The rele- 
vant operators are universal and expressed in terms of 
the light fermion fields and gluons, while their coeffi- 
cients depend on the details of the theory in question. 
A model-independent analysis is complicated by the rel- 
atively large number of contributing operators and by 
the fact that the dominant contributions vary for differ- 
ent models. Furthermore, within a given model, in many 
cases different operators dominate in different regions of 
the parameter space. Heavy paramagnetic systems are an 
exception in this respect: their EDMs are dominated by 
just two terms which are enhanced approximately as Z 3 ; 
one term is directly proportional to the electron EDM 
d e , the other stems from electron-nucleon interactions, 
parametrized by a dimensionless parameter C$- 

In deriving limits for the electron EDM from the corre- 
sponding measurements, commonly firstly the uncertain- 
ties of the numerical proportionality factor are ignored 
and secondly the other relevant term is set to zero, i.e. it 
is assumed that no cancellations occur. When performing 
a quantitative analysis, these issues are obviously impor- 
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tant, especially when keeping in mind that theoretical 
limitations can change the obtained limits by orders of 
magnitude, as observed for the hadronic limits from the 
Mercury (Hg) system, see e.g. [2]. Finally, the obtained 
limits are usually displayed as "allowed" and "forbidden" 
areas in parameter space, making conservative estimates 
obligatory. 



We address both issues in this paper: the first point 
is resolved by more careful estimates for the relevant co- 
efficients below. The second issue can be addressed as 
well, given that at the moment two measurements with 
similar sensitivities are available, from the Thallium (Tl) 
and Ytterbiumfluorid (YbF) systems [3J 0]. However, the 
two systems depend on a similar combination of the two 
terms. Therefore, for the time being, we use in addition 
a limit from a diamagnetic system, namely Mercury [5]. 
While many terms contribute to that EDM, the one ap- 
pearing in paramagnetic systems as well is expected to be 
clearly subdominant; assuming this term to saturate the 
experimental limit therefore constitutes a conservative 
estimate. Together, these three systems allow to obtain 
robust limits for the electron EDM and the coefficient of 
the electron-nucleon interaction, without the assumption 
of vanishing cancellations. 



The outline for this letter is as follows: the second sec- 
tion is devoted to atom EDMs, with a focus on estimates 
of the theoretical uncertainties in their relations to d e 
and Cs- In Sec. Ill, an analogous procedure is carried 
out for molecules, focusing on the EDM of YbF. The 
experimental situation is reviewed in Sec. IV, followed 
by the phenomenological analysis with present data in 
Sec. V, where the new limits on d e and Cs are obtained. 
The results from this analysis allow us to place upper lim- 
its on the EDMs of other paramagnetic systems, which 
we do in Sec. VI, together with an analysis of the future 
prospects, before concluding in Sec. VII. 
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II. EDMS OF ATOMS 

For atoms, Schiff's theorem [5] implies a vanishing 
EDM in the non-relativistic limit for systems of particles 
whose charge distribution is identical to their EDM dis- 
tribution. The limits from the non-observation of these 
EDMs are then related to violations of the conditions for 
this theorem and separated into two classes, depending 
on which of the approximations is more strongly violated. 
For reviews, see e.g. Refs. [7J [5] . 

In paramagnetic atoms, which are our main concern, 
relativistic effects are more important. They are largely 
enhanced for atoms with a large proton number, scaling 
at least like d ~ Z 3 . As mentioned above, this implies a 
sensitivity mainly to the electron EDM, but also a subset 
of electron-nucleon interactions is enhanced. The domi- 
nant component of the latter is described by 

U In=% J2C»(NN)(en 5 a), (1) 

V 1 N=n,p 

where we neglected operators with other Dirac structures 
which are negligible here, but have a largely enhanced 
relative influence in diamagnetic systems. 

In diamagnetic atoms the finite size of the nucleus is 
the main source for the violation of Schiff's theorem. The 
dominant contribution to the corresponding EDM stems 
from its Schiff moment, which is in turn related to quark 
(colour) EDMs and CP-violating four-quark interactions. 
However, the above electron-nucleon interaction is rele- 
vant as well. We will use this fact to obtain an upper 
limit on the coefficient Cs- 

A. The EDM of paramagnetic atoms 

For paramagnetic atoms which have one unpaired elec- 
tron, mainly this electron determines the EDM of the 
atom, as the effects of the ones in closed shells cancel. 
Relativistic effects for atoms with large proton number 
lead to enhancement factors for the electron EDM of 
0(100) in these systems. In addition, the coefficient Cs 
of the electron-nucleon interaction might contribute siz- 
ably. Relating the experimentally observable atom EDM 
to these sources involves complex many-body calcula- 
tions, for which a number of methods exist (for a review, 
see e.g. again [8]), the results of which sometimes span 
large ranges. 

The most sensitive system from that class to date is 
Thallium. Calculations for the enhancement factor yield 
4i/4 ~ [—1041, —179]. One reason for this large range 
is the presence of strong cancellations between differ- 
ent contributions. Recent calculations narrow down this 
range, however, there is some difference remaining [9"rll4|. 
In order to assess its influence, we define two input sets: 
input set I is using the result 4n D —573(20)4 from 
[T5] , which includes the value from . Set II is using the 
value recommended in 10], 4i ^ —466(10)4, which is 



about 20% lower. As both calculations claim a precision 
clearly excluding the other, this situation indicates large 
systematic effects in one or several of the calculations. 

For paramagnetic atoms, the parametrization in 
Eq. ([I]) leads in the limit of infinite nucleon mass to 

%i N = iG F /V2 Z Nl0 l5P N (r) , (2) 

N—n,p 

with the nuclear densities /0./v(r) normalized to unity and 
Z^ denoting the number of the corresponding nucleon 
in the nucleus. Furthermore assuming pn{t) = p(r) and 
abbreviating 1 Cs = J2 N ZmA~ 1 Cq leads to 

U s eN = iG F /V2A <7s7o75PM , (3) 

which is the Hamiltonian typically used in the atomic 
calculations for the corresponding coefficient. It is ob- 
tained in the same kinds of calculations like the one 
for the electron EDM and is plagued by the same can- 
cellations. The most recent results yield 4i(Cs) — 
-7.0(2) x Hr 18 ecmCs [15| and d Ti (C s ) = -4.06(2) x 
10~ 18 e cm Cs [16] . We assign the former to input set I 
and the latter to set II. The combinations with the second 
term read 

4i = -(573 ±20)d e - (7.0 ±0.3) x 10~ 18 ecmC* s ,(4) 
4^ = -(466 ±16)4 -(4.1 ±0.1) x 10~ 18 e cm C* s , (5) 

where we increased slightly the uncertainties of those co- 
efficients which have not been confirmed by an indepen- 
dent calculation. We note that the values for set I are 
consistent with the ratio obtained analytically in [17j . 
while the ratio of the values for set II are significantly 
higher. 

Another interesting system is Cesium (Cs), for which 
several measurements are prepared at the moment, see 
Table [TTJ For this system, the cancellations commented 
upon above are absent, leading to a more stable pre- 
diction. Recent calculations yield compatible results, 
d Cs = (120.5±1.3)d e + (0.801±0.004)xl0- 18 C' s piUE] 
and d Cs = (124 ± 4)4 + (0.759 ± 0.022) x 10" 18 C s ecm 
9J, motivating 

4s = (123 ± 4)4 + (0.78 ± 0.02) x 10~ 18 e cm C s , (6) 

which constitutes in this case an even more conservative 
estimate. This result is consistent with the ratio obtained 
in [13. 

For Rubidium (Rb), the calculations are similarly sta- 
ble, and a very sensitive measurement is prepared as well, 



Note that this definition in principle implies a dependence of Cs 
on the system considered. However, because of (Z n + Z p )/A = 1 
and C'g ~ Cg, this is usually neglected. In addition, the ratios 
Zp,n/A are approximately universal for all atoms considered here 
anyway. 
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see Table [TT| We obtain [TS] 

d Rh = (25.7 ± 0.8)4 + (0.110 ± 0.003) x 10~ 18 ecmC s . 

(7) 

Note that in this case only one recent calculation exists 
for the single coefficients. The uncertainty chosen reflects 
the difference to the analytic ratio given in [T7] and is of 
similar size as the largest difference between experimen- 
tally and theoretically determined CP-conserving quan- 
tities in [18]. 

Finally, there are also plans to measure the EDM of 
the heaviest alkali atom, Francium (Fr), see once more 
Table [il] For this system, even larger enhancement fac- 
tors are expected, d^/de ~ 900 [HI [20]. The coefficient 
of the electron-nucleon contribution has not been calcu- 
lated yet, we use the results of [T7] to estimate its value 
and add an additional 10% uncertainty for that in light of 
the level of agreement for the atoms discussed above. Of 
course a dedicated study of the second coefficient would 
be welcome to confirm this estimate. The result reads 

d Fl = (903 ± 45)4 + (10.9 ± 1.7) x 10~ 18 e cm C s , (8) 

where we conservatively assigned the estimated 5% un- 
certainty in [19] to the coefficient of d e . 



tions are used to obtain the others EH OH [33]. 2 The 
one relevant in this context reads 

C S j o 1.9 x 10 3 (1 + 0.3ZV)" 1 r 2/3 ^ 1 x Ct(<tn) , 

(9) 

where C T {cr N ) = (Cj.{a p ) + C^(cr n )\ , (er p „) implies 

the average over the protons/neutrons in the nuclear 
state and fi denotes the magnetic moment of the nucleus 
in terms of the nuclear magneton /in. We expect the 
uncertainty for these relations to be small, C(%), and 
therefore negligible in this context, as also indicated by 
a recent explicit calculation for a variety of atoms |15j . 
For the tensor coefficient, defined by 

d Hg (C T ) = C^ T x lO-^CV^ecm, (10) 

recent results read C5 S = —5.1 ED and C5 S = —4.3 
[24] . Thus we obtain, using Eq. |9]), 

d Hg (C s ) = -(0.00047 ±0.00005)C S x 10" 18 ecm, (11) 

where we inserted nn g = 0.506 /i n 3 and (crjv) — 
— 1/31/ 1, the estimate from a simple shell model for the 
nucleus, and used the common convention d = dl/I. For 
a more detailed analysis of this system, the reader is re- 
ferred to [2]. 



B. The EDM of Mercury 

For diamagnetic atoms, i.e. atoms with vanishing total 
angular momentum, mainly finite-size effects of the nu- 
cleus determine the EDM. More specifically, its dominant 
source is the CP-odd nuclear Schiff moment [§]. How- 
ever, in the following we will make use of the fact that 
additional sources from clcctron-nuclcon interactions and 
the electron EDM are present. Regarding the latter, 
the value usually used in the literature for Hg reads 
dng(de) = 1-16 x 10~ 2 d e [21]. The corresponding cal- 
culation, however, shows a high sensitivity to higher or- 
der effects; the "corrections" to a previous estimate [2"2"] 
amount to ~ 200% and change the sign. The authors 
point out the sensitivity to correlation effects (which have 
been found to be large for Hg for its other coefficients), 
making a new calculation mandatory. In light of this situ- 
ation we do not see a way to extract a meaningful upper 
limit on the electron EDM from Hg until the theoreti- 
cal situation improves. However, even taking the central 
value quoted above, the bound would be weaker than the 
one from Tl or YbF. 

The electron-nucleon interactions are induced via the 
operators in H e N = 2j=s pt^-In- The coefficients 
in the expression for G?Hg (Cs,p,t) are obtained again 
in atomic calculations; usually only the coefficient of 
the tensor operator is calculated, defined via Hj N = 
G f /V2J2n^t {Ni^a^N^ea^e), and analytic rela- 



III. THE EDM OF PARAMAGNETIC 
MOLECULES 

Polar molecules exhibit very large internal fields, which 
average out to zero in absence of an external field due to 
molecular rotation. The application of an external field 
mixes rotational levels of opposite parity and induces two 
effects: one energy split which is sometimes called some- 
what sloppily an EDM, because it scales as |E ext | for siz- 
able fields, but is T-even, and a much smaller one, which 
is actually T-odd, in which we are interested and which 
is described below. The main difference to atoms is that 
the external field is only used to prohibit the cancellation 
of the effect of the internal field, which is the one acting 
on the electrons. This is why polar molecules can ex- 
hibit huge enhancement factors, increasing the sensitiv- 
ity to fundamental parameters like d e |25j . Analogously 
to atoms, the molecules are categorized according to the 
total angular momentum of their electrons. We discuss 
in the following the paramagnetic case. 

The sensitivity of paramagnetic molecules therefore 
stems in principle from the same mechanism as in para- 
magnetic atoms, but is even higher. As in the case of 



2 Note the different conventions for in different publications, 

,T,P jP.T i \ jT,P ,T,P T/r 

e -g- d atom = d^jtT N) versus d,; nm = dJ^I/I. 

3 Source: WcbElcments (http://www.webelements.com/) 
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atoms, the two main sources are the electron EDM and 
electron-nucleon interactions. Different molecules like 
YbF or PbO are used, which provide a naturally high po- 
larizability. They exhibit effective amplification factors 
of internal versus external fields of O(10 6 ), resulting in 
principle in a sensitivity to the electron EDM of 0(100) 
times that for atoms. 

From the theory point of view, the difficulty lies in 
calculating the relevant internal field, Ej nt , which cannot 
be measured. For this, again multi-body calculations are 
employed, which are complicated by the presence of the 
second core and, as before, the large number of electrons. 
The corresponding interaction energy can be written as 



AE 



lYbF 



\(w dt 



W c Cs) (n-i){E cxt ), 
(12) 

with an external electric field E cxt = E cxt z, n denoting 
the direction of the molecular axis and their alignment 
depending on the external field. The factor 1/2 is due to 
the spin of the electron 4 and the constant Cs has been 
introduced in Eq. In 4J, (n • z)(E eKt ) = 0.558 holds 
|26j . The constant W d j1 reflects the maximal effective 
electric field acting on the valence electron. As noted 
above, in contrast to the atomic case, the effective elec- 
tric field is now given in terms of the internal field, the 
effect of which stops canceling out once the external field 
is applied, due to the closeness of the corresponding ro- 
tational levels. 

Again in parallel to the experimental efforts there has 
been recent theory activity. The relevant results for YbF 
are shown in Table [i As pointed out in [IT], the ra- 
tios of these matrix elements can be estimated analyti- 
cally. Their value for YbF, W d /W c = 114 x 10 18 /ecm, is 
in agreement with the latest calculations [2"Trf2"9"] within 



~ 10% 5 , reflecting the spread in the values for W d - We 
conservatively allow for these 10% variation in both di- 
rections as an error estimate. In absence of a second 
recent determination of W c , we assign it as well as an 
error estimate there, which yields finally 



W d = -(l.l±0.1)xl0 25 Hz/ecm, W c = -(92±9)kHz. 

(13) 

We note that a calculation of W c by a second group with 
the presently available methods would be welcome. From 
these considerations we finally obtain 

d Y bF = -(1.3±0.1)xl0 6 d e -(9740±960)xl0^ 18 ecmC s , 



to be compared with Eqs. Q-Q. 



(14) 



4 Note again the presence of different conventions in the literature: 
Wd is sometimes defined without this factor. 

5 Note that Ref. 1301 is aiming at an analytical estimate rather than 
high precision, they estimate the accuracy to ~ 25%, making 
their result compatible with the following estimate. 



) 25 Hz/ecm) 


W c {kRz) 


Ref. 


Year 


-0.91 


- 82 


EB 


1996 


-1.26 


-120 


SUES 


1994/97 


-1.20 


-104 




1998 


-1.20 


-108 


M 


1998 


-1.21 




m 


1998 


-1.50 




m 


2008 


-1.04 


-92 


nuns 


2007/08 


-1.16 




[29] 


2009 



TABLE I: Calculations for the coefficients in the expression 
for the dipole moment of YbF. 



System Present limit (e cm) Expected limit (e cm) 

199 Hg (0.49 ±1.50) x HT 29 [5] 

205 T1 -(4.0 ±4.3) x 10~ 25 [3] 

133 



Cs* 
Rb 



1.4 x 10" 23 [37] 



1 x 10" 18 [H] 
(1.2 x 10~ 23 ) f [42] 



O(10- 26 /10" 27 ) 38-40 
C(10" 27 /10" 28 ) sm 



210 Fr 

YbF (3.5 ±8.7) x 10" 



e>(Kr 26 /i(T 29 ) [g 

dm] o(io- 22 /io- 23_24 ) 



45 



TABLE II: Present limits on absolute values of EDMs at 
95% CL for the most sensitive atoms/molecules, together with 
short-term/mid-term expected sensitivities. *: Given in the 
paper as (-1.8 ±6.7 ±1.8) x 10 _24 ecm. f : unpublished 



IV. EXPERIMENTAL STATUS 

At present, the most stringent limits relevant to the 
extraction of d e and Cs stem from searches for EDMs of 
Tl g], YbF [J and Hg [5J, see Table [d] Although these 
limits have different orders of magnitude, their different 
dependence on the fundamental parameters of the theory 
actually leads to similar sensitivities. Especially, despite 
the very different factors in Eqs. (|4|/([5]) and (14 1, the 



resulting limits for the electron EDM are similar so far. 
Note, however, that the result for YbF is still statistically 
limited. 

Recently there have been several developments which 
allow to expect significantly improved sensitivities in the 
near future, see also e.g. [51 l46rt48] : the first option is to 
improve the methods described above. With the exper- 
iments for Tl completely dominated by their systematic 
errors, significant advancement seems difficult within this 
system. An improvement, up to two orders of magnitude, 
might come instead from the Cs, Rb and Fr systems [551 - 
|4T)1 14"51 PH] . The expected limits correspond to probing 
the electron EDM to < 10 _29 e cm in the short-term fu- 
ture (1-3 years), and even sensitivities down to 10~ 31 e cm 
seem achievable. 

Further measurements with paramagnetic molecules 
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are expected to strengthen the present limit by another 
order of magnitude or more for YbF, and many more sys- 
tems are explored as well, see e.g. [3E] for a recent list, 
making for an expected improvement of at least the one 
from atoms. 

In the future, trapped molecular ions might also be 
used as sensitive probes for EDMs, however, at the mo- 
ment there are still severe experimental and theoretical 
challenges to overcome. Furthermore, also solid state 
systems are being explored as sensitive probes for the 
electron EDM [33 [SO] ■ While again some experimental 
as well as theoretical progress is necessary before com- 
petitive results can be achieved, recent results show the 
progress in this field |51j . Finally, new techniques are be- 
ing explored for measuring the EDMs of charged particles 
directly by using a storage ring [52 55 . While the main 
focus here is on other systems, there are also proposals 
to use the technique for molecular ions, see e.g. |56j . 

The plethora of ongoing and planned experiments, all 
aiming at the strengthening of present limits by several 
orders of magnitude, will take this field to a new level. 
Especially if one or several of these experiments should 
result in a significant non-zero signal, the question of 
a more refined analysis of the various uncertainties will 
be posed, making a global analysis obligatory. We will 
explore steps in this direction below. 

V. A ROBUST LIMIT ON THE ELECTRON 
EDM 

With the results of the last sections at hand, we 
proceed to derive limits on the electron EDM and the 
electron-nucleon coupling. We do this in two steps: first, 
we derive the limit just from the measurements with 
Tl and YbF, to avoid even input from the conservative 
bound on C$ from Hg. Then we add this as a third con- 
straint, obtaining a much stronger limit on both, d e and 
Cs- We use the data as given in Table [nj i.e. not (yet) 
transforming the given values into symmetric bounds. 

The results are shown in Figs. [I] and [2] for the two input 
sets, where the constraint from each system is shown in 
the d e — Cg-plane. We illustrate by the light grey area in 
Fig. [I] the bound on the electron EDM obtained by the 
combination of the Tl and YbF constraints only (compare 
also to [17]). For input set II, these two constraints are 
basically parallel; therefore no bound can be obtained in 
this case without further assumptions. The dark area in 
the middle is the global fit to all three constraints. The 
projections on the parameters of interest read 

d{ = (0.026 ± 0.065) x 10" 26 ecm and (15) 
& s = (-0.08 ±0.36) x 10~ 7 (16) 

for input set I and 

di 1 = (0.024 ±0.066) x l(T 26 ecm and (17) 
Cg = (-0.12 ± 0.36) x 10~ 7 (18) 




-1.0 -0-5 0.0 05 1.0 

d,/(10- Is e cm) 



FIG. 1: Bounds from Hg, Tl and YbF in the de-Cs-plane, 
using input set I. The very light grey vertical bound indi- 
cates the ID-limit on d e when using only the Tl and YbF 
constraints without the aid of Hg. 
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-1.0 -0.5 0.0 05 10 

d,/(10 _Is e cm) 



FIG. 2: Bounds from Hg, Tl and YbF in the d e -Cs-plane, 
using input set II. 

for input set II. These values are to be compared with 
d e = (-0.31 ± 0.35) x 10~ 26 e cm and C s = (3.2 ± 3.3) x 
10~ 7 , obtained using only the two constraints from Tl 
(input set I) and YbF. The corresponding upper limits 
at 95% C.L. are 

\d{ JI \ < 0.14 x 10~ 26 ecm and (19) 
\& S JI \ < 0.72(0.74) x 10~ 7 (20) 

for the global fit, whereas \d e \ < 0.89 x 10~ 26 ecm and 
\Cs\ < 8.6 x 10~ 7 when excluding the input from Hg. 
Independent of the input set, the global fit therefore re- 
sults in a limit on the electron EDM very similar to the 
one obtained naively from YbF alone, but is obtained in 
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a more conservative manner. Using only paramagnetic 
systems at the moment worsens this limit approximately 
by a factor of six for input set I and is not possible for 
set II. 

In the next section we will show, however, that even 
the conservative assumption entering here via the input 
from the Hg system can be avoided with future data. 



VI. UPPER LIMITS FOR OTHER SYSTEMS 
AND FUTURE PROSPECTS 

The already available limits for the EDMs of Cs and Rb 
given in Table [TT] do not strengthen the limits on d e and 
Cs derived in the last section. This in turn implies that 
we can place non-trivial bounds on these EDMs from our 
results in Eqs.( 15 ),( 16 ). To do so, we map the 95% C.L. 



area from the global fit onto the corresponding interval 
of the atom EDMs, taking additionally the theoretical 
uncertainties there into account. As the obtained inter- 
vals for the two input sets are very similar, we do not 
differentiate in the following between the two. 

Starting with Cs, we obtain with aid of Eq. ^ the 
95% C.L. interval 



d Cs G [-1.6,2.0] x 10" 25 ecm 



(21) 



Therefore the dedicated experiments are expected to im- 
prove the present limit by approximately two orders of 
magnitude before becoming sensitive to possible non- 
vanishing contributions. The same is true for the Rb 
experiments, where the interval reads 



d Rh G [-3.4,4.1] x 10" 2b ecm 



(22) 



Also for Fr we obtain a rather strong limit already, 

d Fl G [-1.3, 1.5] x l(T 24 e cm . (23) 

However, with the expected final sensitivities, see Ta- 
ble [TlJ the planned experiments will be able to improve 
greatly the present bounds or to finally obtain a non-zero 
result. A contradicting measurement in one of these sys- 
tems would indicate a severe issue in one of the involved 
experiments or the theoretical description. 

While the use of the constraint from the Hg system 
proves very advantageous at the moment, in principle it 
would be preferable to perform a similar procedure with- 
out this input. We therefore investigate to what extend 
this is possible with coming measurements, see once more 
Table |Hj To that aim we plot in Figs. [3] and [4] the con- 
straints expected from the future measurements within 
1-3 years and in the longer run, respectively. Note that 
the plotted areas correspond to 1/60 (1/2000) that of 
Fig. [T] The constraints are chosen such that their cen- 
tral value is zero, thereby reflecting the resulting limits 
in the absence of a non-zero result; for significant mea- 
surements, of course all constraints should still overlap. 
Already in the scenario for the first plot the additional 



lH 0-0 




-0.03 0.00 0.05 

d t H10~' s e cm) 



FIG. 3: Bounds from various paramagnetic systems as ex- 
pected in the short-term future (1-3 years) in the d e -Cs- 
plane, see Table [II] Note the different scales compared to 
Figs. [T] and [2] 




-0 005 0.000 
rf,/(10 _i5 e cm) 



FIG. 4: Bounds from various paramagnetic systems as ex- 
pected in the mid-term future (more than 3 years) in the d e - 
Cg-plane, see again Table [n] Note the different scales com- 
pared to Figs. 1 1| >| 



constraint from Hg is rendered unnecessary and the re- 
sults for d e and Cs are significantly improved. These 
plots illustrate clearly the importance of various experi- 
ments with different atoms and/or molecules. First of all, 
at least two competitive measurements are necessary to 
yield a model-independent constraint on d e . Ideally they 
should constrain very different combinations of d e and 
Cs, as for example Rb and Fr. Secondly it is important 
to have more than two constraints in order to confirm 
the theoretical description and safeguard against possi- 
ble systematic issues. Finally, the combination of more 
constraints yields additional precision, which can indi- 
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cate non-vanishing values for d e and Cs earlier. The list 
in Table [TT| indicates that this challenge is met. 

VII. CONCLUSIONS 

Measurements of EDMs are extremely sensitive probes 
of CP-violating phases beyond the SM. They therefore 
have the potential to reveal NP and will continue to 
strongly constrain possible NP scenarios. The experi- 
ments presently planned and constructed will take this 
field to a new level of precision, challenging many models. 
To meet that precision, bounds from these measurements 
should be derived carefully. We have shown in this let- 
ter that it is possible to go beyond the common assump- 
tion of vanishing cancellations, already with present data. 
Doing so, we provided expressions for various systems of 
experimental interest, where we focused on a careful es- 
timate of theory uncertainties. This allowed us to obtain 



more robust limits on the electron EDM and the electron- 
nucleon interaction. Despite the more conservative ex- 
traction, these limits match the more naively extracted 
ones in precision, due to the combination of various mea- 
surements. At the moment the additional input from the 
Hg system is necessary, which is possible with conserva- 
tive assumptions. In the future, even these assumptions 
can be avoided, once strong limits and/or determinations 
of the EDMs from more paramagnetic systems are avail- 
able. 
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